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Gamma ray signatures of ultra high energy cosmic ray 
accelerators: electromagnetic cascade versus synchrotron 
radiation of secondary electrons 



Received: date / Accepted: date 

Abstract We discuss the possibility of observing ultra 
high energy cosmic ray sources in high energy gamma 
rays. Protons propagating away from their accelerators 
produce secondary electrons during interactions with cos- 
mic microwave background photons. These electrons start 
an electromagnetic cascade that results in a broad band 
gamma ray emission. We show that in a magnetized Uni- 
verse (B > 10~ 12 G) such emission is likely to be too 
extended to be detected above the diffuse background. 
A more promising possibility comes from the detection 
of synchrotron photons from the extremely energetic sec- 
ondary electrons. Although this emission is produced in 
a rather extended region of size ~ lOMpc, it is expected 
to be point-like and detectable at GeV energies if the 
intergalactic magnetic field is at the nanogauss level. 

PACS 98.70.Sa ■ 98.70.Rz 

Understanding the origin of the Ultra High Energy 
Cosmic Rays (UHECR) with energy above ~ 5 10 19 eV 
constitute a real challenge for theoretical models, since 
their acceleration requires extreme conditions hardly ful- 
filled by known astrophysical objects. Previous studies 
considered a number of potential sources, including gamma 
ray bursts, active galactic nuclei, large scale jets and neu- 
tron stars, but results are still inconclusive pQ. 

If, as many believe, UHECRs have an extragalactic 
origin, a suppression, known as the Greisen-Zatsepin- 
Kuzmin (GZK) cutoff, should be observed in the spec- 
trum at an energy ~ 7 10 19 eV This is because dur- 
ing propagation energetic protons undergo P7 interac- 
tions (mainly photo-pion production) in the Cosmic Mi- 
crowave Background (CMB) with an energy loss length 
which decreases dramatically from ~ 650 to ~ 20Mpc 
between energies of 7 10 19 to 3 10 20 eV jHJ. The detection 
of a number of events above 7 10 19 eV raised interest in 
non-acceleration models that involve new physics and do 
not require the existence of a cutoff in the spectrum £Q . 
However, the issue of the presence of the GZK feature 
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is still debated |3J. In the following we assume that as- 
trophysical objects capable of accelerating UHECRs do 
exist and we discuss possible ways to identify them. 

If UHECRs are not significantly deflected by the in- 
tergalactic magnetic field (IGMF) their arrival direction 
should point back to the position of their accelerators. 
In principle, this fact can be used to identify accelerators 
located at distances smaller than the proton loss length, 
that constitutes a sort of horizon for CRs. However, even 
in this case, the small event statistics and the high uncer- 
tainty in the determination of the CR arrival direction 
would make any identification problematic. The better 
angular resolution and collection area of AUGER will 
hopefully provide adequate data to address this issue. 

The interactions between UHECRs and CMB pho- 
tons generate secondary gamma rays and electron-positron 
pairs (hereafter referred to as electrons) that in turn ini- 
tiate an electromagnetic (EM) cascade in the univer- 
sal photon background 0. Such a cascade would ap- 
pear to a distant observer as a flux of GeV/TeV pho- 
tons. Recently, Ferrigno et al. [H] suggested to identify 
the sources of UHECRs by searching for this radiation. 
Their calculations show that, in an unmagnetized Uni- 
verse (Bigmf — G), a steady source emitting isotropi- 
cally 2 1 43 erg/s in form of UHECRs (E > 10 19 eV) can 
be detected by a Cherenkov telescope like HESS up to 
a distance of ~ 100 Mpc. In this case the source would 
be point-like, because in the absence of IGMF the EM 
cascade is one-dimensional and propagates radially away 
from the accelerator. Unfortunately, the scenario changes 
dramatically if the more realistic case of a magnetized 
Universe is considered. This is because low energy elec- 
trons produced during the last steps of the cascade are 
effectively deflected and eventually isotropized if their 
Larmor radius is smaller or comparable with the Comp- 
ton cooling length. This condition is satisfied when the 
IGMF is above B iso ~ 10~ 12 (£; 7 /TeU)G, E 1 being the 
energy of the Compton photon. This implies that unless 
the IGMF is extremely weak (<< Bi so ), electrons emit 
Compton photons after being fully isotropized. Thus, an 
extended halo of emitting pairs forms around UHECR 
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sources. The radiation from the halo is emitted isotrop- 
ically, resulting in a very extended, and thus hard to be 
detected, gamma ray source [5]- 

However, if the IGMF close to the accelerator is at the 
nG level, first generation electrons [E ~ 10 19 eV) gener- 
ated during pj interactions cool rapidly by emitting GeV 
synchrotron photons and the development of the cascade 
is strongly inhibited. We propose the possibility to de- 
tect these synchrotron photons from UHECR sources 0. 
This is of great interest because of the following reasons. 
Both synchrotron emitting electrons and parent protons 
are extremely energetic and not appreciably deflected by 
the IGMF, at least on the first Mpcs distance scale. For 
this reason, synchrotron photons are emitted in the same 
direction of parent protons. Thus, they move away from 
the source almost radially, and the observed radiation is 
expected to be point-like, and thus easily detectable and 
distinguishable from the extended cascade component. 
Remarkably, a detection of these sources would allow 
to infer the value of the IGMF close to the accelerator, 
constraining it in the range » 10~ 10 4- 10~ 8 G. Finally, 
since the Universe is transparent to GeV photons, power- 
ful UHECR accelerators located outside the CR horizon 
might be identified in this way. 



which the particle energy is roughly divided in half in 
every collision. This phase ends up with a large num- 
ber of low energy electrons and photons. On the other 
hand, if a IGMF is present, electrons also lose energy 
via synchrotron emission, subtracting energy to the cas- 
cade. In Fig. we show the effective electron loss length 
for Compton/pair production (solid), together with the 
synchrotron loss length for a IGMF equal to 0.1, 1 and 
10 nG (dashed lines). It can be seen that if the IGMF is 
at the level of 1 nG or more, all the electrons with energy 
above ~ 10 18 eV cool fast via synchrotron losses and the 
development of the cascade is strongly suppressed. 

Summarizing, three different regimes, corresponding 
to different values of the IGMF, can be distinguished: 

• Regime I: B IGMF < B lso ~ 10~ 12 G. The EM 
cascade is not affected at all by the IGMF. 

• Regime II: B iso < B IGMF < B syn ~ 10~ 9 G. No 
energy is subtracted to the EM cascade due to syn- 
chrotron losses, but low energy electrons are effec- 
tively isotropized by the IGMF. 

• Regime III: Bigmf > B syn . The development of 
the EM cascade is strongly suppressed since its very 
first steps due to strong synchrotron losses. 



1 Development of the electromagnetic cascade 

In this section we describe the development of the EM 
cascade initiated by a UHECR. Consider a proton with 
energy E P:2 o = E p /10 20 eV. The typical energies of pho- 
tons and electrons produced in pr) interactions are ~ 
10 19 £ p ,20 and ~ 5 10 lf %, 2 o eV respectively 0. In the 
absence of IGMF, such electrons and photons interact via 
Compton and pair production processes with photons in 
the CMB and radio background. In general, this would 
lead to the development of an EM cascade, in which 
the number of electrons and photons increases rapidly. 
In fact, due to the extremely high energy of the parti- 
cles considered here, each interaction occurs in the limit 
r = ebE ^> 1, where e& and E are the energies of the 
background photon and of the energetic electron (pho- 
ton) respectively, both calculated in units of the elec- 
tron rest mass energy. Under this condition the Compton 
scattering is in the extreme Klein-Nishina limit, namely, 
the upscattered photon carries away most of the energy 
of the incoming electron. The same happens during a 
pair production event, in which most of the energy goes 
to one of the two outgoing electrons. Therefore the prob- 
lem reduces essentially to a single-particle problem, in 
which a leading particle loses continuously energy and 
changes state from electron to photon and back due to 
alternate Compton/pair production interactions. Thus, 
the effective loss length of an energetic electron can be 
identified with the loss length of the leading particle 
When the leading particle loses its energy until r sa 1, 
the cascade enters the particle multiplication phase, in 



2 Regime I: one-dimensional cascade 

If the IGMF strength is much less than B iso ~ 10~ 12 G, 
electrons in the cascade do not suffer synchrotron losses, 
nor are they deflected. Thus, the EM cascade develops 
along a straight line. In this case, the calculations by 
Ferrigno et al. apply, and nearby and powerful UHECR 
sources might be detected as point like TeV sources by 
currently operational Cherenkov telescopes [S]. 

In principle, since the strength of the large scale IGMF 
is basically unknown JOjj such a low values of the field 
cannot be ruled out. However, this is probably not a 
good assumption in the vicinity of UHECR accelerators, 
where the IGMF is expected to be appreciable, especially 
if such accelerators are, as it seems reasonable to believe, 
correlated with the structures in the Universe. The cas- 
cade might still be one-dimensional if its last steps de- 
velop sufficiently far away from the source, in a region of 
very low IGMF. Another necessary condition is that the 
IGMF close to the source must be small enough (^; 10~ 9 
G) to avoid a suppression of the cascade due to syn- 
chrotron losses of first generation electrons. 



3 Regime II: extended pair halos 

If the IGMF is strong enough to deflect the electrons 
in the cascade, but not enough to make synchrotron 
losses relevant (namely, 10 _12 G < Bigmf *C 10~ 9 G), 
then the EM cascade fully develops, low energy electrons 
are isotropized, and a very extended pair halo forms 
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around the UHECR source. For an isotropic source, the 
size of the halo can be roughly estimated as follows. 
Let E° s be the energy of the gamma ray photons ob- 
served from the Earth. Such photons are CMB photons 
Compton-upscattered by electrons with energy E e ~ 
20(£° b 7TcV) 1 / 2 TeV. These are the electrons forming 
the pair halo. Since electrons are rapidly isotropized in 
the IGMF, one can assume that they do not propagate 
away from the sites in which they are created. Elec- 
trons in the halo are in turn produced by parent pho- 
tons with energy EP ar < E e . Since the photon mean free 
path against pair production in the infrared background 
X pp decreases rapidly with increasing energy jTH] , we can 
safely neglect the contribution to the halo size from older 
generation (higher energy) photons. Thus, the size of the 
halo can be roughly estimated as lhalo ~ X pp (E^ ir ) 0. 
For a ~ 20 TeV photon the mean free path is about a few 
tens of megaparsecs (see Fig. 2 in JI]). In fact, for the sit- 
uation considered here, the size of the halo is even bigger, 
since the UHECR protons and the first generation elec- 
trons propagate ~ 10 4- 20 Mpc before initiating the EM 
cascade (see Fig. 2). Thus, a conservative estimate of the 
apparent angular size of the halo at 1 TeV can be given 
by: ~ (Ikaio/D) « \Q°(l halo /2QMpc){D /imMpc)" 1 , 
where D is the distance of the source. This indicates 
that pair halos are extremely extended, bigger than the 
field of view of Cherenkov telescopes (the HESS field 
of view is ~ 5°) and thus hardly detectable. The problem 
becomes even worse if one considers also protons with en- 
ergy below 10 20 eV, which have much longer loss length, 
increasing up to ~ 1 Gpc for proton energies equal to 
~ 5 10 19 eV. On the other hand, such protons are likely 
to contribute only to the TeV flux of very distant sources. 

In the recent work by Armengaud et al. [H] the deflec- 
tion of electrons in the IGMF has been neglected, even 
when a IGMF stronger than Bi SO ~ 10 -12 G was as- 
sumed. On the other hand, the authors considered the 
deflection of ~ 10 20 eV protons, which is in fact totally 
negligible if compared with the full isotropization of elec- 
trons. Therefore their claim about the detectability of 
cascade gamma-rays seems to us over-optimistic. 

The cascade emission peaks at TeV energies (SJ, mak- 
ing a detection by GLAST problematic. 



4 Regime III: synchrotron gamma rays 

If the IGMF close to the UHECR accelerator is at the 
level of 1 nG or above, the development of the cascade 
is strongly suppressed, since the very high energy elec- 
trons produced during p"f interactions cool rapidly via 
synchrotron losses before undergoing Compton scatter- 
ing. It is evident from Fig. 2 that for a ~ nG IGMF this 
is true for electron energies well in excess of E e ~ 10 18 
eV. Such electrons emit synchrotron photons with en- 
ergy E syn w 2( J B/nG)(£/10 19 eV) 2 GeV, detectable by 
GLAST. It is important to stress that our results are sen- 
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Fig. 1 Effective electron loss length. Solid: Compton/pair- 
production losses in the CMB and radio background (with a 
cutoff at 2MHz 12;). Dashed: synchrotron losses. 



sitive only to the value of the IGMF close to the source, 
while they are unaffected by the value of the field on 
much larger scales. This is because synchrotron emit- 
ting electrons are produced within a proton interaction 
length l pl w 10 Mpc from the accelerator. As a conse- 
quence, the only assumption required is that the size of 
the magnetized region surrounding the accelerator must 
be greater or comparable with ^ P7 .Superclusters of galax- 
ies constitute an example of large and magnetized regions 
satisfying our requirement |10j . 

We now estimate the angular size of the synchrotron 
emission. After propagating over an interaction length, 
a proton of energy E p is deflected by an angle d p « 
0.8°(10 20 eV/£ p ) (B /nG) v /(/ P7 /10Mpc) ^/(Zc/Mpc) where 
l c is the IGMF coherence length ^3]. Due to the high 
energies considered, secondary electrons produced in jry 
interactions move in the same direction of the parent 
protons. In a cooling time electrons are deflected by an 
angle t? e ~ aX/RL, Rl being the electron Larmor ra- 
dius and a a number of order unity representing the 
probability that the leading particle is actually an elec- 
tron [S]. Remarkably, if expressed as a function of the 
synchrotron photon energy, the deflection angle is inde- 
pendent on the magnetic field strength and reads: i? e ~ 
0.5 o (E syn /10GeV)~ 1 . Thus, an observer at a distance D 

would see a source with angular size: i) b s ks \j^p + (% 

that, for D = 100 Mpc and for photon energies of 1 4 10 
GeV is of the order of a fraction of a degree. This is 
comparable with the angular resolution of GLAST, that 
would classify these sources as point-like if they are lo- 
cated at a distance of ~ 100 Mpc or more. This leads to 
the important conclusion that, even if synchrotron pho- 
tons are produced in an extended region of size ~ l pl sur- 
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Fig. 2 Spectra for a source located at lOOMpc. The luminos- 
ity in UHECRs is 2 10 44 erg/s, with spectral index is S — 2. 
TOP: B IG mf = 0.5(curve 1), 5 (2), 50nG(3), E cut = 10 21 eV. 
BOTTOM: E cut = 5 10 20 , 10 21 , 510 21 eV, B IG mf = lnG. 
Dotted: intrinsic spectra. Solid: spectra after absorption. 

rounding the accelerator, the resulting gamma ray source 
would appear point-like to a distant observer. 

If the proton spectrum extends well above 10 20 eV, 
these sources might also be detected by Imaging Atmo- 
spheric Cherenkov Telescope arrays (IACT) operating at 
energies above 100 GeV. The angular resolution of these 
instruments is a few arcminutes, and thus the sources will 
appear extended. However, it has been proven that IACT 
are powerful instruments to image extended gamma ray 
sources Jl] , and thus they might still detect and map the 
emission from UHECR sources. Finally, powerful UHECR 
accelerators located at a distance of 1 Gpc or more would 
appear as point sources. In the next section we discuss 
the energy requirement for a detection. 



5 Detectability and energetics 

Fig. |3 shows synchrotron spectra for a steady source at 
a distance of 100 Mpc in a uniformly magnetized re- 
gion of size 20 Mpc. Since the radiation is produced 
in a region about one interaction length away from the 
source, the actual structure of the field is not crucial. 
Steady state proton and electron spectra have been cal- 
culated taking into account all the relevant energy losses 
and proton escape from the magnetized region. Solid 
lines have been computed taking into account the opac- 
ity of the Universe to very high energy photons due to 
pair production in the cosmic infrared background (TBI, 
while dotted lines show the unabsorbed spectra. The to- 
tal luminosity in UHECRs with energy above 10 19 eV is 
Luhe = 2 1 44 erg/s, with a differential energy distribu- 



tion cx E~ s exp(—E/E cut ), with 8 = 2. Results are quite 
insensitive to the slope of the spectrum. 

In the top panel of Fig. |2] E cut = 10 21 eV and the 
IGMF is equal to 0.5, 5 and 50 nG (curves 1, 2 and 3). If 
the IGMF is significantly greater than ~ 50 nG, the peak 
of the emission falls at TeV energies, where absorption is 
very strong. On the other hand, if the field is well below 
~ 0.5 nG, synchrotron emission becomes unimportant 
and the cascade contribution dominates. However, for 
the broad interval of values of the IGMF strength between 
0.5 and 50 nG ; the formation of a synchrotron point-like 
gamma ray source seems to be unavoidable. 

In the bottom panel of Fig. our predictions are 
compared with the sensitivities of GLAST and of a generic 
IACT such as HESS or VERITAS. A IGMF of 1 nG is 
assumed and the different curves refer to values of the 
cutoff energy in the proton spectrum equal to 5 10 20 , 10 21 
and 5 10 21 eV (top to bottom). For such a field, the con- 
dition for the detectability of a point source by GLAST 
is roughly L UHE > 8 10 43 -. 2 1 44 (D/100Mpc) 2 erg/s for 
5 = 2.0 -.2.6. In contrast, for IACT arrays the minimum 
detectable luminosity is roughly 2 orders of magnitude 
higher, since the source has to be located at a distance 
of ~ lGpc in order to appear point-like. However, less 
powerful accelerators can still be detected as extended 
sources. Since the peak of the emission falls at ~ 10 GeV, 
future IACT operating in the energy range 10-. 100 GeV 
would be powerful tools to search for these sources. 

If the CR spectrum smoothly extends down to GeV 
energies with slope 8 — 2, the required total CR lumi- 
nosity for a source to be detected by GLAST is Lcr > 
5 1 44 (D/100Mpc) 2 erg/s. For a beamed source, the re- 
quired luminosity is reduced by a factor /& ~ 0.02($b/10°), 
i?6 being the beaming angle, and the detectability con- 
dition reads: L CR > 10 43 (/ 6 /0.02)(£>/100Mpc) 2 erg/s. 
This luminosity is small if compared, for example, with 
the power of an AGN jet, that can be as high as 10 47 erg/s 
|lfij . Thus, astrophysical objects that can in principle 
satisfy the energy requirement for a detection do exist. 
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